Do Topological Charge Solitons Participate in DNA Activity? 
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We present a novel electromagnetic model of DNA molecules, in which the P-bonds act as tunnel 
junctions and the H-bonds as capacitors. Excess charge in the model gives rise to two coupled mod- 
ified sine-Gordon equations, which admit topological solitonic excitations. We study the dynamics 
of the solitons, their effect on the DNA transport properties, and comment about their role in the 
DNA functioning. We propose specific experiments in order to test our predictions. 
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The DNA molecule is central to every living organ- 
ism. It carries the genetic code of the organism, and 
continuously controls the synthesis of proteins, which are 
vital to its functioning and existence. Since the pioneer- 
ing discoveries of Watson and Crick, a considerable ef- 
fort was devoted to the study of the DNA, and much 
has been understood about its structure and activity |Q . 
Yet there are some key features which pose open ques- 
tions. One such feature is the long range correlation and 
control between segments of the DNA. Synthesis of pro- 
teins is done via local "reading" of specific sectors in 
the DNA. The initiation of the reading is done by an- 
other segment, which can be located far away along the 
DNA sequence (thousands of bases away) . Consider this 
and other examples Q , it seems as if the DNA molecule 
has the capability of transmitting information over long 
distances and in a specific manner (the information is 
transmitted to/from specific targets). Possible efficient 
candidates for such a transmission are solitons and soli- 
tary waves (rather than point particles/wave packets and 
linear waves). The existence and propagation of confor- 
mational solitons related to the DNA replication process 
have been studied in the past H . Non-topological charge 
solitons in proteins have also been studied Q . Motivated 
by our studies of topological charge solitons in one dimen- 
sional (ID) arrays of mesoscopic tunnel junctions ||, ||, 
we have investigated the possibility of the existence of 
topological charge solitons in DNA. 

We have developed a novel electromagnetical model of 
the DNA molecule, which is based on the properties of 
charge dynamics in the DNA (see Fig. |l|). We view each 
of the DNA strands as a ID array. A unit (or 'grain') 
of the array is a sugar and a base attached to it. The 
grains are connected longitudinally by Phosphate groups 
(P-bonds), and parallely, i.e., to the other strand, by the 
Hydrogen bonds (H-bonds) between bases. An additional 
electron residing on an atom belonging to the sugar-base 
grain can hop from atom to atom, thus obtaining a ki- 
netic energy. The H- and P-bonds form barriers to the 
charge propagation. The proton in the H-bond can ef- 
fectively screen a net charge density on either side of the 
bond by shifting its position towards this side. (By 'net' 



charge we mean the change from the charge distribution 
of the unperturbed DNA.) As a result, electrons do not 
cross the H-bond. Hence the bond can be viewed as a 
capacitor. The P-bond barrier stems from the two oxy- 
gens which are transversely connected to the phosphorus. 
These oxygens share three electrons with the phosphorus, 
giving rise to two a bonds and one tt bond. As the tt elec- 
tron can be shared with both oxygens, it behaves as an 
electron in a double well, and occupies the lowest level. 
When another electron approaches the well it encounters 
a barrier due to the energy difference to the next level 
of the well. However, since this barrier is narrow, the 
approaching electron can tunnel through the well. Thus 
from the charge dynamics point of view, the P-bond be- 
haves as a tunnel junction. 

We start with the model for a single strand (see Fig. |l|). 
For simplicity we assume that all the bases are of the 
same type. Each grain i is composed of four sub-grains 



(j = 0, 1, 2, 3), assigned with phase variables, 



These 



phases are related to the electric potential of the sub- 
grains through 4>ij (t) = i J_ Vij(t')dt'. The conjugate 
variable to each </>i.j is the charge on the sub-grain, Qij. 
The Lagrangian of this model is 



C 



k i 12, ^0 \2 , 



2L/2 



'H,l - <Pi,2 



2L/2 



'H,2 - <Pi,3 ) 



EjCOs(0i+i,i - (pi,3) 



L /2 



*H,0 - 9i,2 



(1) 



It includes three types of energies: inductive energies 
which represent the hopping of electrons, capacitive en- 
ergies which represent the capacitive properties of the H- 
and P-bonds, and a tunneling energy, which represents 
the tunneling process in the P-bond. This latter energy 
is proportional to the cosine of the phase difference across 
the bond, according to the tight binding picture JtJ . The 
values of the parameters can be obtained, in principle, 
from experiments, but for now we employ only a quali- 
tative view. The Lq/2 inductance denotes the hopping 
from the open side of a base to its sugar, while the L/2 
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inductance denotes the hopping between P-bonds. Since 
the former involves more hopping around carbon rings, 
we assume that L > L. Both Cs and Co denote in this 
single-strand model capacitances to the outside world. 
Cs is the capacitance seen from the sugar, while Co is 
the capacitance seen from the open H-bond. We assume 
that both are of the same order of magnitude and much 
smaller than C, the capacitance of the P-bond tunnel 
junction. Ej is the tunneling strength of the P-bond. 
The length scale in our model, a, is the distance between 
grains, which is 3.4A. (In the model we use units in which 
a = 1.) 
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FIG. 1. A schematic image of a DNA molecule (left), and 
the double-strand model (right). On the left, P denote the 
P-bonds between the sugars S and the dashed lines are the 
H-bonds between the bases A,G,T and C. The parameters 
of the double-strand model are described in the text. In the 
single-strand case only one half of the model is used, and Co 
denotes the capacitance to the outside world. 

From the kinetic part of (|l|) we see that out of the four 
variables describing the i'th grain only three are inde- 
pendent. According to our experience in the study of ID 
arrays of tunnel junctions ||, we introduce non-local (or 
integral) charge variables. These variables are very use- 
ful in trying to probe the collective, non-local dynamics 
of the DNA. We envision a propagation of charge from 
a certain point in the far left of the chain ('— oo') to a 
certain point in the far right ('oo'), passing on its way the 
i'th grain. Referring to the i'th junction as the junction 
between the i'th and the (i + l)'th grains, we define Ri 
as the charge that has passed through this junction: 



Ri = J]] Qi'J — + Qi' ,2 + Qi>,0 + Qi') ~ Qoo 

i+1 

= Qi,i + 5^(Qf, 2 + Qi',o) ■ (2) 
■i+i 

Si+i i Qi' ,j means summation over the sub-grain 
charges, starting from (i' = i + l,j = 1) and going right, 
with the internal order (i+1, 1) — > 2) — > (i+1, 0) — > 

(i + 1, 3) — > (i + 2, 1), and so on. Next we define as the 
charge that has reached junction i: 



Qi'J = Qi,9 + & = Z)(Qi',a + Qi\o) ■ (3) 



After the canonical transformation to the variables: 7T; = 
4>i+i,2 - <f>i,2 - Qi, 0i = 0-i+i,i - 0i,3 and 4> ifi = fafi - 
4>i,2, we obtain the following Euclidean Lagrangian (in 
the continuum limit which we justify later): 
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where Cso = CsCo/(Co + Cs) is the effective capacitance 
across the chain (parallel coupling of Cs and Co). As 
Cs is much smaller than Co, Cso ~ Cs. We define a 
tunneling inductance Lj = ^/(2ttEj) 2 ($ = h/e is 
the flux quantum), which we take to be much smaller 
than L. From (|j) one can identify three characteristic 
frequencies in the system: ft 2 = l/(LC) of the q mode, 
ujj = \/{LjC) of the R mode, and lo 2 = 2/(L C S0 ) 
of the qo mode. In the limits we are working, we have 
£! 2 <C uj 2 ,oj 2 . Thus we can integrate out the R and go 
modes, and obtain an effective theory for the q mode. 
The effective real time Lagrangian is 
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where Ec = e 2 /(2C). The equation of motion of 
the following modified sine-Gordon equation: 
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is 
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where Vd = e/(2irC). This is a Kirchoff law for the 
equivalent electrical circuit. Since the characteristic 
length scale in a sine-Gordon model is C/(Co + Cs), 
the continuum limit is justified when C > Co + Cs, 
which indeed corresponds to our assumption. Transform- 
ing into dimensionless space, time and charge variables: 
x — > x 1 = yJC/Cox, t — > t' = \/LCt, q — * q' = j^q, we 
obtain: 
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where C = L f§ L , and ^ s = C s /C . 

Next we introduce the double-strand model (see again 
Fig. [l]). We distinguish between the two strands by the 
superscripts a,f3. The capacitance Co describes now the 
capacitive coupling between the two strands, i.e., over 
the H-bonds. We assume that it is much larger than the 
capacitance to the outside world, Cs, but still smaller 
than C. We take the parameters of the two strands to 
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be equal. Following the same steps of derivation as in 
the single-strand model (details will be given elsewhere 
H ) , we obtain two (dimensionless) coupled modified sine- 
Gordon equations of motion: 
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rfe = , (9) 



where \i = Cs/2C*o, and we have used the average and 
the relative charge variables: 



(10) 



Similar equations, without the £ term, where studied in 
the past in connection with the stacked Josephson junc- 
tions model H . 

We turn now to study the modified sine-Gordon equa- 
tions we have derived, starting from the single-strand 
model. The pure Sine-Gordon equation has exact topo- 
logical solitons solutions. We have checked numerically 
and found that the extra term in equation (^) does not 
effect the stability of the topological soliton solution, 
though it does induce interaction with the plasmons. 
This interaction causes the soliton to slowly radiate away 
its kinetic energy. Using the collective coordinate: 



X = / q x xdx, 

Zn 



(11) 



we find that the soliton velocity actually oscillates with 
relatively high frequency, but its average is almost con- 
stant, decreasing only slightly over large period of time. 
Therefore, though the soliton's motion is not persistent 
in the exact sense, charge can propagate relatively large 
distances along the strand without the need to apply a 
driving force in the form of a potential difference. 

The extra term in equation (Q) destroys the Lorenz 
invariance of this equation. The dispersion relation for 
small amplitude linear waves is 
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and the group velocity tend to zero for both small and 
large k values. When the soliton velocity exceed the max- 
imal group velocity, it leaves in its wake all the small 
amplitude waves. 

In the double-strand model, the dimensionless param- 
eter yU serves as a measure of the strength of the interac- 
tion between the strands (it is small for strong interac- 
tion). Since C$ <C Co, fi is very small. In the extreme 
fj, = limit, q x should be zero as well, in order that the 
Lagrangian for this model would be finite. As both q a 
and qP have integer values at the two edges of the chain, 



q « 2im everywhere. We are thus left with a single equa- 
tion for q: 
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This equation is equivalent to the one of the single-strand 
model, but now q represents a soliton anti-soliton pair 
(with zero total charge) rather than a single charged soli- 
ton. 

Equations (||,^|) have two simple topological solutions. 
One is the "symmetrical" solution, namely, q = and q 
is a usual Sine-Gordon kink. This solution is unstable ,at 
least at low energy, since it corresponds to two charges 
with equal sign which tend to separate. It is, however, 
possible that it becomes stable at high energy, as was 
demonstrated for the stacked Josephson junctions model 
fllPf . The other "anti-symmetrical" solution has q = 
and q in the form of a kink. It is a stable solution which 
correspond to an electron-hole pair (each on a different 
strand). It is plausible that this excitation is used in the 
DNA to transmit information over long distances. It can 
be created at a specific segment, which is responsible for 
the initiation of a certain protein synthesis, by tunneling 
of an electron through the H-bond (as a result of exter- 
nal agent). The pair can, then, propagate almost freely 
along the DNA, and be annihilated at the protein synthe- 
sis segment by a reverse tunneling, thus transmitting the 
order for this synthesis. The targeting of a specific an- 
nihilation segment is probably done through the specific 
base sequence of both segments, which was neglected in 
our model. 




FIG. 2. Topological solution representing stationary charge 
in the a strand. The values /1 = 0.1 and £ = 1 where used. 
The steepest curve is q(x), and the other is q(x). Numerical 
results are indicated by continues lines and the theoretical 
curves of equation (|l4|), by dashed lines. 

There is also a way to approximate a solution which 
corresponds to an electron (hole) moving with velocity v 
in the, say, a strand: 



q = 2 arctan exp [—j{x — vt) / c] 
q = 2 arctan exp [—7(2; — vt) / c] 



(14) 
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where 7 is the relativistic factor, and c = 1/ \/2(fi + 1), 
5=1/ V^Mi are the linear wave velocities in equations 
(|H|) (when the last term in equation (J9J) is neglected). 
They correspond to the different linear waves velocity in 
the two equations. The exact solution can be obtained 
numerically and is shown in Fig. |^. 

An important way to test our model is via measure- 
ments of the current- voltage characteristics. An electron 
can be injected to one of the strands to form a topolog- 
ical charge soliton, and may be subjected to an applied 
external voltage V and an Ohmic dissipation R (both per 
unit length) . We implement this in our model by adding 
the terms F + aq to equation (||) and aq to equation 
@, where F = V/Vd and a = RC/y/LC, are the exter- 
nal force and the dissipation rate respectively. For given 
values of force and dissipation rate a soliton reaches a 
limiting velocity which corresponds to the measured cur- 
rent. The I-V curves which we obtain for different choices 
of the /1 and £ parameters are shown in Fig. ||. It should 
be noted that above some threshold voltage close to the 
maximum voltage indicated in each curve of Fig. [| the 
system becomes unstable, as more and more soliton anti- 
soliton pairs are created. 
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FIG. 3. The numerical current (velocity) - voltage (force) 
curves for a charge soliton propagating in one strand. In the 
lower curves a — 0.5 and [i = 0.5, while in the upper ones 
a = 0.5 and /1 = 0.05. In the continuous lines ( — 0.1, and in 
the short dash lines £ = 1. The theoretical curve of equation 
( p"5| ) is indicated by a long dash line. 

We proceed to estimate the limiting velocity. To do 
so, we equate the energy gain by the external force and 
the dissipative energy [[[1J. We find that X (7/0 + 7/0) = 
irF/2a. This reduces in the small /1 limit to 
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This theoretical curve is also shown in Fig. |[ The 
ascending of the numerical curve in the "relativistic" 
regime above the theoretical curve is obviously due to 
the £ term. 



We turn now to discuss some possible experimental 
tests of our predictions. The most straightforward test 
is to measure the I-V characteristics of a DNA that was 
mentioned earlier. Such experiments are now being de- 
veloped |jl2). We propose to use first a single strand vs 
a double strand. It is easier to inject solitons into the 
former. In the case of a double strand, it would be useful 
if each end is made of a short single strand (a differ- 
ent one on each end). In both cases it would be easier 
to compare the experimental results with the theoretical 
predictions if artificial strands composed of repetitia of 
one base are used. Another approach would be to study 
the magnetic response of circular strands. The idea is to 
place many circular strands on a surface, apply a time- 
dependent magnetic flux (say, of a saw tooth form), and 
measure the response of the system. For short circular 
strands at low temperature we expect to observe persis- 
tent current carried by the charge solitons. Finally, a less 
trivial experiment, is to inject a charge soliton at the end 
of the DNA via Scanning Tunneling Microscope (STM) , 
and look for a response at the other end (for example, 
via an attached molecule that has a fluorescent response 
to an incoming charge). 

These are just three examples of many more possible 
tests of our predictions about the existence of topolog- 
ical charge solitons in DNA molecules. If turned to be 
verified, we expect charge solitons to have crucial role in 
the DNA activities as means of transfer of information 
and energy over long distances, and to specific locations. 
To study this role, one can use our model as a start- 
ing point, incorporating into it the inhomogeneity of the 
chain due to the different bases, as well as interactions 
with external molecules at specific sites. 

We are most thankful to E. Braun for sharing with 
us his ideas about experimental measurements of DNA 
electrical transport. This research is supported in part 
by a GIF grant. 
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